Benchmark calculation of n- 3 H and p- 3 He scattering 
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The n- 3 H and p- 3 He elastic phase-shifts below the trinucleon disintegration thresholds are cal- 
culated by solving the 4-nucleon problem with three different realistic nucleon-nucleon interactions 
(the I-N3LO model by Entem and Machleidt, the Argonne Vis potential model, and a low-fc model 
derived from the CD-Bonn potential). Three different methods - Alt, Grassberger and Sandhas, Hy- 
perspherical Harmonics, and Faddeev-Yakubovsky - have been used and their respective results are 
compared. For both n- 3 H and p- 3 He we observe a rather good agreement between the three different 
theoretical methods. We also compare the theoretical predictions with the available experimental 
data, confirming the large underprediction of the p- 3 He analyzing power. 



PACS numbers: 21.45.+V, 21.30.-x, 24.70. +s, 25.10.+S 
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I. INTRODUCTION 

The four-nucleon (4N) system has been object of in- 
tense studies in recent years. In first place, this system 
is particularly interesting "theoretical laboratory" 

to test the accuracy of our present knowledge of the 
nucleon-nucleon (NN) and three nucleon (3N) interac- 
tions. In particular, the effects of the NN P-waves and of 
the 3N force are believed to be larger than in the A = 2 or 
3 systems. Moreover, it is the simplest system where the 
3N interaction in channels of total isospin T = 3/2 can 
be studied. In second place, there is a number of reac- 
tions involving four nuclcons which are of extreme impor- 
tance for astrophysics, energy production, and studies of 
fundamental symmetries. As an example, reactions like 



He + 7 or p 



5 Hc 



(the hep 



process) play important roles in solar models and in the 
theory of big-bang nucleosynthesis. 

Nowadays, the 4N bound-state problem can be nu- 
merically solved with good accuracy. For example, in 
Ref. [l[ the binding energies and otherproperties of the 
a-particle were studied using the AV8' 21 NN interaction; 
several different techniques produced results in very close 
agreement with each other (at the level of less than 1%). 
More recently, the same agreement has also been ob- 
tained considering different realistic NN+3N force mod- 
els . 

In recent years, there has also been a rapid advance in 
solving the 4N scattering problem with realistic Hamil- 
tonians. Accurate calculations of four-body scattering 
observables have been achieved in the framework of the 
Alt-Grassberger-Sandhas (AGS) equations @-[H||, solved 
in momentum space, where the long-range Coulomb in- 



teraction is treated using the screening and renormal- 
ization method [H, [H] . Also solutions of the Faddeev- 
Yakubovsky (FY) equations in configuration space [13 - 
[IH and the application of the Hyperspherical Harmonics 
(HH) expansion method 11011 to the solution of this prob- 
lem have been reported |2fj| . l2lj ] . 

In addition to these methods, the solution of the 4N 
scattering problem has been obtained also by using the 
resonating group model (RGM) method [22H25j . Calcu- 
lations of scattering observables using the Green's func- 
tion Monte Carlo method are also underway [26[. 

The 4N studies performed so far have evidenced several 
discrepancies between theoretical predictions and exper- 
imental data. Let us consider first n- 3 H elastic scatter- 
ing. Calculations based on NN interaction models dis- 
agree 0, 0, H3| rather sizeably with the measured to- 
tal cross section [23], both at zero energy and in the 
"peak" region (E n ss 3.5 MeV). This observable is found 
to be very sensitive to the NN interaction model 0]. 
At low energy, the discrepancy is removed by includ- 
ing a 3N force fixed to reproduce the triton binding en- 
ergy 0, HO, [H, [28[ , but it remains in the peak region. 
The analysis of the differential cross section has shown 
similar discrepancies, but definitive conclusions arc diffi- 
cult to extract since the experimental errors are rather 
large. 

In this respect, the p- 3 He elastic scattering is more in- 
teresting since there exist several accurate measurements 
of both the unpolarized cross section [29M3TI ] and the 
proton analyzing power A y o [31H331 ] . The calculations 
performed so far (with a variety of NN and NN+3N in- 
teractions) have shown a large discrepancy between the- 
ory and experiment for A y0 [1, [13, EH, [3J, H3] • In ad- 
dition, at the Triangle Universities Nuclear Laboratory 
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(TUNL), there has been recently a new set of accurate 
measurements of other j>- 3 He observables (the 3 He ana- 
lyzing power Aq v and some spin correlation observables 
as A yy , A xx , A xz , A zx , and A zz ) at E p = 1.60, 2.25, 4.05 
and 5.54 MeV, which has allowed a phase-shift analysis 
(PSA) [35| . A preliminar y co mparison with these data 



has been reported in Ref. [21 



In order to have definite answer about the ability of the 
different interaction models to reproduce the experimen- 
tal data it is certainly of interest to establish the accuracy 
reached by the theoretical methods in the solution of the 
A = 4 scattering problem. In a previous benchmark, the 
results obtained by different groups working with differ- 
ent techniques were found to be at variance with each 
other [l?}. Clearly, this situation should be clarified be- 
fore questioning the ability of present NN+3N force mod- 
els to describe the experimental data beyond the binding 
energy of 4 Hc. This is the purpose of the present paper, 
in which we present low energy n- 3 H and p- 3 He scatter- 
ing results obtained by three different groups, using in- 
dependent methods to solve the four-body problem, i.e., 
the AGS equations, the variational HH expansion, and 
the FY equations. 

The potentials used in this paper are the I-N3LO model 
by Entem and Machleidt HU, with cutoff A = 500 MeV, 
the Argonne V\s (AV18) potential model [37]], and a 
low-fc model derived from the CD-Bonn potential f38j . 
The I-N3LO potential has been derived using an effec- 
tive field theory approach and the chiral perturbation 
theory up to next-to-next-to-next-to-leading order. The 
AV18 potential is a phcnomenological potential having a 
rather strong repulsion at short intcrparticle distances. 
The low— fc potentials have been obtained separating the 
Hilbert space into low and high momentum regions and 
using the renormalization group method [38j to integrate 
out the high-momentum components above a cutoff A. 
The low— fc potential adopted in this work is obtained 
starting from the realistic CD-Bonn potential [3!| and 
using a smooth cutoff A — 2.5 fm _1 . The cut of the 
high-momentum part is reflected in configuration space 
in an almost total absence of the repulsion at short in- 
tcrparticle distances. Note that the first and third model 
are non-local, while AV18 is local in configuration space. 
The three potentials reproduce equally well the np and 
pp data, and are a representative set of the large variety 
of modern NN potential models. We note finally that I- 
N3LO and AV18 interactions, without the inclusion of a 
suitable 3N interaction model, largely underestimate the 
4 He binding energy B( 4 Hc). On the contrary, with the 
adopted low-fc potential model we have B( 4 He) = 29.04 
MeV, slightly overestimating the experimental value of 
28.30 MeV. 

This paper is organized as follows. In Section [Til a 
brief description of the methods used for this benchmark 
is reported. In Section IIII1 a comparison between the 
results obtained within the different schemes is shown. 
In Section ITV1 the theoretical calculations are compared 
with the available experimental data. The conclusions 



will be given in Section [V] 



II. METHODS 

In order to solve the 4N scattering problem we em- 
ploy the AGS equations, the HH method, and the FY 
equations. The various procedures are briefly described 
below. 

The total kinetic energy, T c m ., in the center of mass 
(cm.) and the nucleon kinetic energy, En (N = p, n), 
in the laboratory reference frame are given by 



2n ' 



E 



N 



(1) 



where p, = (3/4)A/jv is the reduced mass of the 1 + 3 
system, Mjv is the nucleon mass, and q the magnitude of 
the relative momentum between the two clusters. 



A. AGS Equations 

The AGS equations [4(| are integral equations for the 
four-body transition operators. They are well-defined 
only with short-range potentials. Nevertheless, together 
with the screening and renormalization method (8|, l4lj , 
they can be applied also to the reactions involving 
charged particles. In the 4N system we use the isospin 
formalism and solve the symmetrized form of the AGS 
equations 0- In this case there are only two distinct 
four-particle partitions, one of the 3 + 1 type and one of 
the 2 + 2 type. We choose those partitions to be (12,3)4 
and (12) (34) and denote them in the following by a = 1 
and 2, respectively. The corresponding transition op- 
erators Up a for the initial states of the 3 + 1 type, as 
appropriate for the n- 3 H and p- 3 He scattering, obey the 
integral equations 

Un = —(GoTGo) 1 P34 — -P34 Ui Go TGo U\\ 

+U 2 G Q TG Q U 21 , (2) 

U21 = (Go TGo)" 1 (I-P34) 

+(1-P 34 )U 1 G TG U 11 . (3) 

Here Go = (E + ie — Hq)^ 1 is the free resolvent, E being 
the energy of the 4N system and Ho the free Hamiltonian, 
and Pij is the permutation operator of particles i and 
j. The (short-range) two- nucleon potential V s enters the 
AGS equations via the two-nucleon transition matrix T = 
V s + V s GqT and the 3+1 and 2+2 subsystem transition 
operators 



U a — P a G 1 + P a TGo U a , 



(4) 



where Pi = P 12 P23 + P13 P23 and P 2 = P13 P 2i - The 3+1 
elastic scattering amplitudes are given by (p/|7"|pi) = 
3(^1 (p/) |^ii l^i (Pi )} where the factor 3 results from the 
symmctrization and \^> a (pj)) arc properly normalized 
initial/final channel state Faddccv components. 
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In order to include the Coulomb interaction V c be- 
tween the protons in the p- 3 He scatterin g w e use the 
screening and renormalization approach pi l4l| . We add 
to the nuclear pp potential the screened Coulomb one 
V R (r) = V c (r) exp(-(r/i?)™). Thus, the AGS equa- 
tions with V s + V R are well-defined but all transition 
operators and the resulting amplitudes depend on the 
screening radius R. The renormalization procedure in 
the R -+ oo limit yields the full p- 3 He transition ampli- 
tude 



(p/l^c)|p,) = (P/l^ m 1p i > 



lim (p f \[T {R) -t c R m -}\ Pi )Z- f 



(5) 



where (p/|t^ m |pi) and (pf\t% m '\pi) are the proper and 
screened Coulomb amplitudes between the cm. of two 
charged clusters, respectively; the former is known an- 
alytically. The renormalization factor is defined in 
Ref. Q. Thus, the long- and Coulomb-distorted short- 
range parts in the scattering amplitudes are isolated and 
their infinite R limit is calculated separately. The long- 
range part of the amplitude (p/|£ij m '|pi) is of two-body 
nature and its R — > oo limit after renormalization is 
(pf\t^j m -\pi). The Coulomb-distorted short-range part 
[7(h) — t^ m -] is calculated by solving the AGS equations 
for V s + V R numerically at a finite R that is sufficiently 
large to get i?-indepcndcnt results after the renormal- 
ization. In other words, the R -+ oo limit is reached 
with sufficient accuracy at finite R. However, R must 
be considerably larger than the range of the nuclear in- 
teraction thereby leading to a slower partial-wave con- 
vergence. The right choice of the screening, i.e., the 
exponent n, is essential in dealing with this difficulty. 
For a fast convergence with R we have to ensure that 
V R (r) approximates well V c (r) for r < R and simulta- 
neously vanishes smoothly but rapidly for r > R, provid- 
ing a comparatively fast convergence of the partial- wave 
expansion. Using the optimal value n = 4 we obtain 
reasonably converged results with R ranging from 10 to 
15 fm and including two-proton partial waves with or- 
bital angular momentum up to 10. The R- convergence is 
slower at lower energies; the worst cases are the S waves 
at E p = 2.25 McV where we estimate the accuracy of 
our phase-shift results to be around 1%. In contrast, the 
?i- 3 H results are converged very well, considerably better 
than 0.2%, as demonstrated in Ref. Q where also the 
details on the included partial waves can be found. 



B. HH Expansion 

The wave function describing a n- 3 H or p- 3 He scatter- 
ing state with total angular momentum quantum num- 
bers J, J z , incoming relative orbital angular momentum 
L, and channel spin S (S = 0,1) can be written as 



LS,JJ Z 
1+3 



LS,JJ Z 
C 



LS,JJ Z 



(6) 



where the part ' JJz describes the system in the re- 
gion where the particles are close to each other and their 
mutual interactions are strong. Hence, \$jLS,JJ z VSin [ s \y CS 
in the limit of large inter-cluster distances. This part 
of the wave function is written as a linear expansion 
12 (i c jl SJ y(L i wnere is a set of basis functions con- 
structed in terms of the HH functions (for more details, 
see, for example, Ref. fl9l]). 

LS J J 

The other part \\! A ' z describes the relative motion of 
the two clusters in the asymptotic regions, where the 1+3 
interaction is negligible (except eventually for the long- 
range Coulomb interaction). In the asymptotic region 

LS J J LS J J 

the wave functions \l/ 1+ 3 z reduces to ^ A ' % which 
must therefore be the appropriate asymptotic solution of 
the Schrodinger equation. Let us consider, for example, 
the p- 3 He case. Then, ^/ A ' JJ * can be decomposed as a 
linear combination of the following functions 



lc LS,JJ z 



JJ Z 



^2[Yl(vi)® ® s lls 

i=i 



qyi 



qyi 



where yi is the distance between the proton (particle I) 
and 3 Hc (particles ijk), q is the magnitude of the relative 
momentum between the two clusters, si the spin state of 
particle I, and </>3 is the 3 He wave function. Moreover, Fl 
and Gl are the regular and irregular Coulomb function, 
respectively, with 77 = 2^,e 2 /q. The function ft{y) = 
[1 — exp(— (3y)] 2L+1 in Eq. has been introduced to 
regularize Gl at small y, and /i(y) — > 1 as y is large, 

LS J J 

thus not affecting the asymptotic behavior of \& 1+ 3 z . 
Note that for large values of qyi , 



fL{yi)G L {r),qyi) ±iF L (r),qyi) -> 

exp ±i(qyi - Ln/2 - rjh\{2qy{) + a L ) 



(8) 



where ol is the Coulomb phase-shift. Therefore, Q~^s jj 
(^ls J J, ) describe the asymptotic outgoing (ingoing) 
p- 3 Hc relative motion. Finally, 



LS.JJ. 



E 

L'S' 



Sll'S S s'^ ls ,jj z 



qJtt q+ 
°LS,L'S' i! 'L'S',JJ z 



(9) 

where the parameters Sjjg L , s , are the S-matrix elements 
which determine phase-shifts and (for coupled channels) 
mixing parameters at the energy T c m .. Of course, the 
sum over V and S' is over all values compatible with 
the given J and parity n. In particular, the sum over L' 
is limited to include either even or odd values such that 
(-1) L ' = (-) L = n. 

The S-matrix elements S^g L , s , and coefficients 

LS J J 

occurring in the HH expansion of ^ c ' * are determined 
by making the functional 



[S. 



LS,L'S' 



o 



LS.L'S' 



E\ * 



LS,JJ Z 
1+3 



(10) 
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stationary with respect to variations in the S[g L , s , and 
c ls J (j£ onn variational principle) . In the above equation, 
E = r c . m . — £?( 3 Hc) is the energy of the system, B( 3 Hc) 
being the 3 He binding energy. By applying this principle, 
a linear set of equations is obtained for S'[g L , s , and c£ . 
This linear system is solved using the Lanczos algorithm. 

This method can be applied in either coordinate or 
momentum space, and using either local or non-local po- 
tentials [l9| (see also Ref. [42[ for an application to the 
A = 3 scattering problem). The first step is a partial 
wave decomposition of the asymptotic functions f2 jj , 
a task which can be rather time consuming, in particular 
for the J 11 = 2~ state. After this decomposition, the cal- 
culation of the matrix element in Eq. (fTU| is fast. Then, 
the problem reduces to the solution of the linear system, 
which is performed using an iterative method (however, 
this solution has to be repeated several times due to the 
necessity to extrapolate the results, see below). 

The expansion of the scattering wave function in terms 
of the HH basis is in principle infinite, therefore a trun- 
cation scheme is necessary. The HH functions are es- 
sentially characterized by the orbital angular momentum 
quantum numbers ii, i = 1,3, associated with the three 
Jacobi vectors, and the grand angular quantum number 
K (each HH function is a polynomial of degree K). The 
basis is truncated to include states with ^1+^2+^3 < ^ m ax 
(with all possible re-coupling between angular and spin 
states appropriate to the given J). Between these states, 
we retain only the HH functions with K < K maK . In 
the calculation we have included only states with total 
isospin T = 1. 

The numerical uncertainty comes from the numerical 
integrations needed to compute the matrix elements of 
the Hamiltonian and the truncation of the basis. It has 
been checked that the numerical uncertainty of the cal- 
culated phase-shifts related to the numerical integration 
is small (around 0.1 %). The NN interaction has been 
limited to act on two-body states with total angular mo- 
mentum j < j max = 8 (at the considered energies, greater 
values of j max are completely unnecessary). The largest 
uncertainty is thus related to the use of a finite basis. 
The convergence with £ max is rather fast and the value 
4iax = 6 have been found to be sufficient. The main 
problem is related to the slow convergence of the re- 
sults with K max . This problem can be partly overcome 
by performing calculations for increasing values of K max 
and then using some extrapolation rule (see for example 
Ref. HH) to get the "-ftT max — > 00" result. This procedure 
has an uncertainty which can be estimated. A detailed 
study of this problem will be published elsewhere [43| . 
The convergence of the quantities of interest in term of 
-Kmax is slower when NN potentials with a strong repul- 
sion at short interparticle distance are used such as for 
the AVI 8 potential. In this case we have estimated the 
uncertainty to be of the order of 0.5 % in the extrapolated 
phase-shifts. This problem is less relevant for the I-N3LO 
and the low-fc models. In these case, the uncertainty has 
been estimated to be at most 0.3 %. 



C. FY Equations in Configuration Space 

In late sixties, Yakubovsky [44| has managed to gener- 
alize the three-body equations derived by Faddeev [45( to 
an arbitrary number of particles. These equations were 
primarily derived for a system of particles submitted to 
short range pair- wise potential V s . Nevertheless it be- 
comes possible to include also repulsive Coulomb inter- 
action if these, from now on called Faddeev- Yakubovsky 
equations, are formulated in configuration space. To this 
aim, we split the Coulomb potential V c into two parts 
(short and long range), V c = V S ' C + V . The splitting 
procedure is quite arbitrary, one should only take care 
that the long range part V LC of the Coulomb potential 
approaches sufficiently fast the full Coulomb interaction 
V c when any of interparticle distances becomes large. 
The simplest application of FY equations is the prob- 
lem of four identical particles. They result into a set of 
two differential equations coupling the two so called FY 
components, namely Kf 2 3 and and have the form: 



E-Ho-VZ-^Vtf 



(V{ 2 + V{f)Pt [(1 



-"■12.3 



^34)^2,3 



E 



l.G 

13 



(V12 + V 1 S 2 C )P 2 [(1 



rr34 
H 12 



^34)^2,3 



,(H) 



,(12) 



where P± , P 2 and P34 are the particle permutation oper- 
ators, equivalent to those described in Section Hi A[ and 
e = ±1 is a phase accounting for the Pauli principle be- 
tween two identical particles (e = +1 for bosons and 
e = — 1 for fcrmions). 

Each FY component F = (K, H) is considered as a 
function of its proper set of Jacobi [r| vectors x,y,z 
and expanded in angular variables for each coordinate 
according to 



x 



xyz 



The quantities F a are called regularized FY amplitudes 
and Y a are tripolar harmonics, containing spin, isospin 
and angular momentum variables. The label a holds for 
the set of 10 intermediate quantum numbers describing 
a J x , T = 1 state in the partial wave basis. 

The FY components F — (K 1 H) are subject to 
Dirichlet-type boundary condition imposed on a three 
dimensional rectangular grid. Both components vanish 
on any of three (x, y, z), well as at the borders 

X = X max and y — Umax 

of the chosen resolution grid. 
On contrary, a boundary condition equivalent to Eq. (|9]) 
is imposed on the z — z max border for the FY compo- 
nents of type K; if no 2+2 particle channels are open FY 
component of type H must also vanish at the z — z max . 
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As discussed in the previous section, the expansion of 
the scattering wave function in terms of the partial wave 
basis is in principle infinite and a truncation scheme is 
necessary. In this work the partial wave basis was trun- 
cated to include all the states with j x < 4, j„ < 4 and 
j z < 3, in the so-called j — j coupling scheme |5L"ll5l]. By 
studying the convergence of the calculated phase-shifts 
with respect to the size of the partial wave basis, we 
have concluded that this truncation scheme should pro- 
vide results accurate at 1% level. 

The numerical implementation of these equations is 
described in detail in Ref. 1151. 



III. RESULTS 

In this section we present the phase-shifts for the most 
relevant waves calculated using the three different meth- 
ods described above. The selected energies for n- 3 H arc 
E n = 1, 2, 3.5 and 6 MeV, while for p- 3 He are E p = 2.25, 
4.05 and 5.54 MeV, corresponding to cases where exper- 
iments have been carried out . 

The states considered are those with J v = ± , 1 ± 1 
and 2~. The scattering in other J v states is dominated 
by the centrifugal barrier and therefore the phase-shifts 
arc smaller and not very sensitive to the interaction and 
the method used to calculate them. Note that, for the 
J n = 2~ state, we have chosen to report only the 3 P-2 
phase-shift, since the 3 F 2 phase-shift and the relative 
mixing parameter are in any case very small. Never- 
theless the coupling between the 3 P2 and 3 F 2 waves has 
been included in the calculations, since the presence of 
the 3 i 7 2 component in the asymptotic part of the wave 
function has a sizable effect on the 3 P 2 phase-shift. 

Let us remember that the S-matrix for elastic n- 3 H 
and p- 3 He scattering has dimension 1 for = ± states 
and dimension 2 for J > 0. In the first case, the S'- 
matrix is parametrized as usual as S[g LS = exp(2i<5^g). 
For J > 0, since the S'-matrix is unitary and symmetric, 
we can write it as 



S = o T s n o 



with Sd a diagonal matrix written as 
{Sd)ls,l'S' = Sll'Sss'£ 2iSl 



(14) 



(15) 



where 8jjg is the phase-shift (in the Blatt-Biederharn rep 
resentation) of the wave LS. Due to the unitarity prop 
erties, 8jjg is a real number. The matrix O in Eq. |T 
is parametrized as 



O = 



cose 



Jk 
c Jtt 



sine 



Jtt 
.Jtt 



where e Jv is the so called mixing parameter of the given 
state. Clearly the values of the phase-shifts and mix- 
ing parameters may depend on the (arbitrary) choice on 
the coupling scheme between the spin of the two clus- 
ters and the spherical harmonic function Yl(d) in the 



asymptotic functions flf^ s (see, for example, Eq. ([7])). It 
can be shown that the phase-shifts defined as discussed 
above are independent on such choices, while the mixing 
parameter, on the contrary, depends on them. Neverthe- 
less, it is easy to establish the linear relation to trans- 
form the mixing parameter from one coupling scheme to 
another. In the following, we chose to report the mixing 
parameters defined in the LS coupling scheme by Eq. ([7]). 
Moreover, in the following tables, the values reported in 
a column labeled as 2S+l Lj (using the "spectroscopic 
notation") are relative to the phase-shift S'[^. 

In Table U we present the phase-shifts, mixing parame- 
ters, and total cross sections for n- 3 H scattering obtained 
using the I-N3LO potential at the selected energies. By 
inspecting the table, we can notice the good agreement 
between the three different techniques. The maximal de- 
viation of the results is less than 1%, fully in line with the 
estimated errors. Furthermore, the agreement between 
the results of AGS and HH techniques is even better, 
only in a few cases the HH and AGS results differ by more 
than 0.5%. The strongest deviation, of the order of 0.4 
deg, is observed with the FY results at the largest stud- 
ied energy. This slightly larger deviation might be due to 
the necessity by the FY method to perform the transfor- 
mation of the aforementioned potential to configuration 
space. In this respect, we note that the AGS calcula- 
tion is fully performed in momentum space, while in the 
HH calculation, part of the needed matrix elements are 
calculated in momentum space (those involving &q ' JJ * , 
which have to be calculated with more accuracy), and 
part in configuration space (those involving ^^ S ' JJ '). 

In Table [II] we present the same n- 3 H results obtained 
using the AV18 potential. In this case, the convergence 
of the HH expansion is more problematic, in particular 
due to the necessity to extrapolate the HH results. We 
observe that we still have a very good agreement for the 
1 So, 3 Po, 3 P2, and 1 + phase-shifts, while the differences 
in the l~ phase-shifts appear to be more enhanced. 

In Table HIT1 we have reported the phase-shifts obtained 
with the low— k potential derived from CD-Bonn. In this 
case, the calculations has been performed using the AGS 
and HH methods, only. We again observe an overall good 
agreement between the results obtained by the two tech- 
niques, except for the lowest energy where the differences 
are sizeable. 

The total cross sections at are found to be in agreement 
within 0.05 b at all considered energies. By comparing 
the values obtained using the different potentials, we can 
observe the following well-known characteristics: (i) at 
low energies the I-N3LO and AV18 models overpredicts 
the experimental cross section. For example, at E n = 1 



( 16 ) MeV, a t expt « 1.6 b |23, while a]^ 



r AV18 



1.8 

b. On the contrary, with the low-fc potential the cal- 
culated <r t is quite close to the experimental one. This 
behavior is related to the strict relation between the to- 
tal cross section at low energy and the triton binding 
energy [3, [H [H [H| . (it) at the peak (around E n = 3.5 
MeV), the experimental cross section has been measured 
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to be (T( Xpt » 2.45 b [27| . In this case we note that the 
AV18 and low-fc potential models underpredict sizeably 
the experimental value, while crJ _N3LO is quite close to 
it. 

Let us now consider p- 3 He scattering. The phase-shifts 
and mixing parameters obtained within the three meth- 
ods have been reported in Tables IIV1 fVl and IVI1 corre- 
sponding respectively to the I-N3LO, AV18, and the low- 
fc NN potential models. Here the differences between the 
various techniques are larger than in the n- 3 H case, espe- 
cially at low energy and for the ,V = 0* states. For the 
AV18 potential, we note that the HH results arc morcless 
intermediate between the AGS and FY results. 

In Table [VI] we have also reported the phase-shifts and 
mixing parameters obtained by the recent PSA [35| . Note 
that the low-fc potential used in this work is the only 
potential which docs not underestimate the 4 He bind- 
ing energy. The PSA estimates have rather large errors. 
However, it is possible to draw some conclusions about 
the capability of this low-fc potential model to describe 
the experimental data. As can be seen, the PSA S-wave 
phase-shifts seem to be well reproduced (except for the 
1 5'o phase-shift at 5.54 MeV) by the calculations. Also 
the 3 Po and 1 P\ agree well, but for these cases the exper- 
imental errors are large. On the other hand, we note a 
sizeable underestimation of the large 3 Pi and 3 P2 phase- 
shifts. 

Let us now see how the fairly good agreement found for 
the phase-shifts and mixing parameters calculated with 
the three different methods reflects on the observablcs. 
We have considered the differential cross section and the 
neutron (proton) analyzing power A v q for n- 3 H (p- 3 He) 
elastic scattering at the considered energies, as functions 
of the cm. scattering angle. Furthermore, we have also 
considered the triton ( 3 He) analyzing power A$ y . This 
observable is in fact rather sensitive to small variations 
of the phase-shifts in the kinematical regime considered 
in this paper. 

In Figs. [T] and [2] we have reported the results obtained 
using the AGS equation (solid lines), the HH expansion 
method (dashed lines), and the FY equations (dotted 
lines) using the I-N3LO potential. As can be seen by 
inspecting the two figures, the three curves almost al- 
ways perfectly coincide and cannot be distinguished. We 
have also reported the experimental data for the n- 3 H 
differential cross section [46| and the three p- 3 He observ- 
ables [291 - I331 . [35| . We note that the differences between 
the three calculations, where they can be appreciated, are 
in any case always smaller than the experimental errors. 

The agreement between the three calculations when 
the AV18 potential is adopted is again rather satisfac- 
tory, as can be seen in Figs. [3] and [U A small disagree- 
ment can be observed only for the Aq v observable (see 
the panels in the last row of FigfJJ. This observable 
is also rather sensitive to the small D-wave and F-wave 
phase-shifts not reported in Tables HT1 and fVl We already 
know that the AV18 model contains a stronger repulsion 
at short interparticle distance than the I-N3LO. As dis- 



cussed above, the convergence of the HH method for this 
case is more problematic and consequently the calcula- 
tions have a larger uncertainty. In spite of these difficul- 
ties, the agreement in the considered observables is still 
quite good. 

Let us consider now the low-fc potential, which has no 
repulsion at short interparticle distance. Consequently, 
in this case, we expect a good agreement between the 
results of the different techniques. For this potential, the 
calculations have been performed using the AGS (solid 
curves) and HH (dashed curves) methods, only, and the 
corresponding results arc reported in Figs. [5] and |H] The 
two curves are practically indistinguishable, confirming 
that for soft potentials the convergence of the calculations 
is excellent. 

Finally, in the literature for p- 3 He scattering, there ex- 
ist measurements of other spin correlation observablcs 
(Ayy, A xx , A zz , A xz , and A zx ). Also for these observ- 
ables we have found a good agreement between the pre- 
dictions obtained by the three different methods, for all 
the potential models considered here. The comparison of 
the theoretical predictions and the experimental data for 
these observablcs will be discussed in the next section. 



IV. COMPARISON WITH EXPERIMENTAL 
DATA 

In this section we discuss the comparison between the 
theoretical calculations and the experimental data. We 
consider here only p- 3 He scattering since for this process 
the experimental data are more abundant and precise. 
The figures presented in this section can be considered 
as an update of previous comparisons 0, H, H3, HH H3 ■ 
For the three observablcs considered so far (da/dCl, A y o, 
and Aq v ), the comparison between theory and experi- 
ment can be inferred already from Figs. [THU However, 
in order to better appreciate the differences in the predic- 
tions obtained by the three potential models as compared 
to the experimental data, we summarize again in Fig. [7J 
the results for da/dfl, A y o, and Aq v . In order to take into 
account the (slight) different predictions obtained using 
the three different theoretical methods, we have decided 
to present the calculated observables for each potential as 
bands. Each band contains the results obtained by using 
the three different methods. As can be seen from Fig. [71 
the differential unpolarized cross sections obtained us- 
ing the I-N3LO potential (red bands) agree well with the 
experimental data. With the other two potentials we ob- 
serve some disagreement, in particular around # c . m . « 30 
deg and in the large scattering angle region. The results 
obtained for A y o are found to depend on the potential 
model. Here, we observe the well known underprediction 
of the experimental data by the theoretical calculations. 
Interestingly, the results obtained with the low-fc poten- 
tial are in a better agreement with the experimental A v q. 
A similar situation is found also for Aq v , as can be seen 
in the three lower panels of Fig. [JJ It is worthy to note 



7 



that the effect of supplemen ting; the AV18 potential with 
the Urbana 3N force model |47| has been found to be al- 
most negligible for this observable [3l[. The inclusion of 
the new chiral 3N potential derived in Ref. [48| is under 
study [lil (see Ref. [2(| for a preliminary report). 

In Fig. HI we report the results found for the A yy and 
A xx spin correlations at the three different proton en- 
ergies. As can be seen, for these two obscrvables the 
predictions obtained with the three potentials are almost 
identical. We observe that the calculated A yy is slightly 
at variance with respect to the experimental data, while 
the A xx observable is reasonably well reproduced by the 
calculations. 

Finally, in Fig. |9]we compare the results obtained for 
the A xz , A zx , and A zz spin correlation observables. In 
this case, only the E p = 5.54 proton laboratory energy is 
considered, since only for this energy experimental data 
exist. Also in this case, the sensitivity to the different 
potential models is not significant. Moreover, the calcu- 
lations reproduce well the (few) experimental data. 

V. CONCLUSIONS 

In this work, we have studied several low energy 
n- 3 H and p- 3 He elastic observables by using three dif- 
ferent approaches, the HH, AGS and FY techniques. 
Around four years ago, some of the authors of the 
present paper presented very accurate solutions of the 4- 
nucleon scattering problem using the AGS technique 0- 
Q. They were able to take into account the long-range 
Coulomb interaction using the screening-rcnormalization 
method [HI, HH. In recent years, also the accuracy of the 
calculations performed using the HH and FY techniques 
increased 0, US HH- Therefore, it becomes appropriate 
to compare the results obtained by the different methods 
in order to test their capability to solve the 4N scattering 
problem. This is the primary aim of the present paper. 
Another important motivation is to provide a set of solid 
converged results in the literature, which could represent 
useful benchmarks for future applications in A = 4 scat- 
tering. 

In the present paper we have shown that for I-N3LO 
and the selected low-fc potential model, which have a 
"soft" repulsion at short intcrparticle distances (the low- 



k model has no repulsion at all), the results obtained 
by the different techniques are in very good agreement. 
With the AV18 potential, the agreement is not so perfect, 
although the (slight) differences can be appreciated only 
for some small polarization observables. We can conclude 
therefore, that the A = 4 scattering problem is nowadays 
solved with a very good accuracy, better than 1%. 

Concerning the comparison with the experimental 
data, we have confirmed the large undcrprcdiction of the 
p- 3 Hc A y o observable, a problem already put in evidence 
some time ago [l4], HH, HI| , and certainly related to the 
N — d il A y puzzle" . For this observable we have observed 
a moderate dependence on the considered potential mod- 
els. The low-fc potential is found to give a better descrip- 
tion of the observable when compared with the experi- 
mental data. However, the same potential does not re- 
produce well the unpolarized cross section. We have also 
found a small undcrprcdiction of the theoretical Aq v and 
A yy observables, while other measured observables, such 
as A xx , A XZl A zx , and A zz spin correlation coefficients, 
show less sensitivity to the potential models. They are in 
good agreement with the available (sparse) experimental 
data. 

The discrepancies found, in particular for A v q, indicate 
a serious difficulty of the existing NN force models in 
describing the 4N continuum. This difficulty can hardly 
be solved by the inclusion of a standard type 3NF, used 
to reproduce the few-nucleon binding energies (l7l . I2H 
l3lj . Its origin could rather lie either in the NN forces 
themselves, or in the presence of a 3NF of unknown type. 
Clearly, an eventual solution of the A = 4 A y o problem 
should be related in some way to the solution of the N — d 
11 A y puzzle" . 

Finally, we conclude noticing that it would be inter- 
esting to extend the present analysis to p- 3 H, n- 3 He and 
d-d scattering observables, which have already been cal- 
culated in the framework of the AGS equations for dif- 
ferent NN interactions 0, ■ 
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En 


'S 


3 Po 


3 S 1 


3 A 


e 1+ 


'Pi 


3 Pi 


e 1 - 


3 P 2 






1.0 


-38.10 
-38.02 
-38.31 


4.15 
4.10 
4.00 


-33.32 
-33.31 
-33.56 


-0.09 
-0.08 
-0.11 


-0.23 
-0.22 
-0.24 


5.99 
5.86 
6.13 


9.63 
9.64 
10.13 


9.44 
9.14 
9.6 


8.98 
8.95 
9.16 


1.77 
1.77 
1.81 


AGS 
HH 
FY 


2.0 


-51.93 
-51.98 
-52.34 


10.54 
10.50 
10.54 


-45.66 
-45.72 
-45.99 


-0.36 
-0.35 
-0.39 


-0.44 
-0.43 
-0.50 


13.13 
13.12 
13.55 


24.18 
24.25 
25.15 


9.15 
9.18 
9.62 


23.96 
23.96 
24.52 


2.13 
2.13 
2.19 


AGS 
HH 
FY 


3.5 


-65.54 
-65.66 
-66.15 


20.31 
20.26 
20.62 


-57.99 
-58.08 
-58.40 


-0.91 
-0.91 
-0.91 


-0.72 
-0.72 
-0.79 


20.74 
20.94 
21.17 


40.94 
40.97 
41.50 


9.45 
9.55 
9.33 


43.98 
43.91 
44.42 


2.38 
2.38 
2.41 


AGS 
HH 

FY 


6.0 


-80.53 
-80.57 
-80.98 


32.71 
32.55 
33.40 


-71.75 
-71.79 
-71.93 


-1.77 
-1.80 
-1.81 


-1.16 
-1.15 
-1.22 


26.88 
26.92 
27.05 


52.35 
52.25 
52.00 


10.62 
10.68 
10.71 


60.04 
60.01 
59.96 


1.97 
1.97 
1.97 


AGS 
HH 
FY 


TABLE I: n- 3 H phase-shifts and mixing 
at E n = 1.0, 2.0, 3.5, and 6.0 MeV. 


; parameters (in degrees) and total cross section a t (in 


barns) for the I-N3LO potential 


E n 


'S 


3 Po 


S S 1 




e 1+ 


'Pi 


3 Pi 


e'~ 


3 P 2 






1.0 


-38.52 
-38.44 
-38.55 


4.36 
4.26 
4.36 


-33.67 
-33.57 
-33.75 


-0.10 
-0.09 
-0.09 


-0.24 
-0.21 
-0.28 


6.15 
5.87 
6.14 


9.64 
9.44 
9.62 


9.38 
9.19 
9.45 


8.94 
8.82 
8.93 


1.80 
1.78 
1.81 


AGS 
HH 
FY 


2.0 


-52.43 
-52.41 
-52.55 


10.93 
10.82 
10.92 


-46.08 
-46.04 
-46.23 


-0.38 
-0.37 
-0.37 


-0.46 
-0.42 
-0.47 


13.30 
13.00 
13.36 


23.90 
23.39 
23.86 


8.99 
9.19 
9.07 


23.45 
23.21 
23.44 


2.12 
2.10 
2.13 


AGS 
HH 
FY 


3.5 


-66.12 
-66.14 
-66.23 


20.75 
20.61 
20.62 


-58.48 
-58.53 
-58.66 


-0.93 
-0.95 
-0.94 


-0.75 
-0.72 
-0.77 


20.73 
20.68 
20.75 


40.09 
39.63 
39.98 


9.24 
9.48 
9.31 


42.51 
42.22 
42.37 


2.33 
2.32 
2.33 


AGS 
HH 
FY 


6.0 


-81.03 
-81.05 
-80.95 


32.77 
32.61 
32.53 


-72.19 
-72.40 
-72.22 


-1.78 
-1.87 
-1.86 


-1.22 
-1.20 
-1.24 


26.53 
26.55 
26.58 


51.13 
51.27 
50.95 


10.37 
10.57 
10.47 


57.87 
57.94 
57.57 


1.93 
1.93 
1.92 


AGS 
HH 
FY 



TABLE II: Same as Table |U but for the AV18 potential. 
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'S 


3 Po 


3 Si 


3 A 


e 1+ 


'Pi 


3 Pi 


e 1 - 


3 P 2 






1.0 


-36.39 


3.52 


-32.03 


-0.08 


-0.19 


5.34 


8.86 


9.79 


8.37 


1.62 


AGS 




-36.08 


3.41 


-31.88 


-0.06 


-0.19 


5.01 


8.70 


9.34 


8.23 


1.60 


HH 


2.0 


-49.73 


9.03 


-43.99 


-0.32 


-0.37 


12.05 


22.61 


9.76 


22.79 


1.96 


AGS 




-49.61 


8.94 


-43.95 


-0.28 


-0.37 


11.83 


22.51 


9.69 


22.71 


1.95 


HH 


3.5 


-62.94 


17.75 


-56.01 


-0.82 


-0.63 


19.72 


39.30 


10.24 


43.20 


2.26 


AGS 




-63.06 


17.74 


-56.10 


-0.79 


-0.63 


19.88 


39.41 


10.32 


43.25 


2.27 


HH 


6.0 


-77.57 


29.44 


-69.51 


-1.66 


-1.03 


26.38 


51.44 


11.57 


60.41 


1.94 


AGS 




-77.77 


29.46 


-69.64 


-1.70 


-1.04 


26.56 


51.48 


11.63 


60.45 


1.94 


HH 



TABLE III: Same as Table [H but for the low-momentum potential derived from the CD Bonn potential. In this case, only the 
AGS and HH results are reported. 





'So 


3 Po 


3 a 


3 ^1 


e 1+ 


'Pi 


3 Pi 


e 1 - 


3 P 2 




2.25 


-40.64 
-41.23 
-41.57 


8.04 
7.73 
7.74 


-35.00 
-35.47 
-35.49 


-0.24 
-0.34 
-0.28 


-0.53 
-0.54 
-0.58 


10.64 
10.42 
10.84 


17.29 
17.11 
17.75 


8.61 
8.69 
8.43 


16.26 
16.11 
16.41 


AGS 
HH 
FY 


4.05 


-58.23 
-58.61 
-59.12 


17.94 
17.76 
18.12 


-50.79 
-51.01 
-51.15 


-0.94 
-0.97 
-0.96 


-0.84 
-0.82 
-0.94 


18.90 
18.97 
19.26 


35.50 
35.43 
35.78 


8.73 
8.85 
8.62 


36.61 
36.53 
36.88 


AGS 
HH 
FY 


5.54 


-68.28 
-68.50 
-69.00 


25.41 
25.07 
25.81 


-60.02 
-60.11 
-60.03 


-1.45 
-1.51 
-1.40 


-1.08 
-1.07 
-1.18 


23.05 
23.00 
23.16 


44.54 
44.34 
44.13 


9.28 
9.36 
9.28 


48.53 
48.29 
48.33 


AGS 
HH 
FY 


TABLE IV 


: p- 3 He ph 


ase-shifts and 


mixing 


; parameters (in 


degrees 


;) for the I-N3LO potential at E v 


= 2.25, 4.05, and 5.54 MeV. 



E p 


'So 


3 Po 


3 Si 


3 Di 


e 1+ 


'Pi 


3 Pi 


6 1 - 


3 P2 




2.25 


-41.11 


8.46 


-35.26 


-0.26 


-0.56 


10.93 


17.35 


8.43 


16.29 


AGS 




-41.53 


7.84 


-35.65 


-0.34 


-0.48 


10.33 


16.76 


8.56 


15.76 


HH 




-41.70 


7.82 


-36.01 


-0.29 


-0.52 


10.52 


17.08 


9.01 


15.90 


FY 


4.05 


-58.70 


18.41 


-51.17 


-0.93 


-0.87 


18.89 


34.83 


8.46 


35.65 


AGS 




-58.93 


17.89 


-51.34 


-0.98 


-0.82 


18.85 


33.49 


8.79 


35.33 


HH 




-59.02 


17.61 


-51.76 


-1.00 


-0.82 


18.57 


34.81 


8.82 


35.36 


FY 


5.54 


-68.75 


25.82 


-60.41 


-1.43 


-1.12 


22.91 


43.65 


9.00 


47.09 


AGS 




-68.96 


25.05 


-60.78 


-1.55 


-1.10 


22.89 


43.20 


9.32 


46.72 


HH 




-68.92 


24.74 


-60.91 


-1.55 


-1.06 


21.93 


44.01 


8.08 


46.53 


FY 



TABLE V: Same as Table 1171 but for the AV18 potential. 
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Ep 


'So 


3 Po 


3 Si 


3 Di 


e 1+ 


'Pi 


3 Pi 


e 1 - 


3 P 2 




2.25 


-38.74 


6.85 


-33.60 


-0.21 


-0.46 


9.63 


15.95 


9.13 


15.13 


AGS 




-38.96 


6.53 


-33.80 


-0.26 


-0.43 


9.21 


15.76 


8.98 


14.91 


HH 




-39.1±1.7 


5 ± 6 


-34.5±0.7 






8±2 


17±4 


10±20 


16.5±0.7 


PSA 


4.05 


-55.74 


15.60 


-48.93 


-0.87 


-0.75 


17.79 


34.64 


9.43 


35.19 


AGS 




-56.09 


15.48 


-49.15 


-0.85 


-0.74 


18.02 


33.93 


9.54 


35.33 


HH 




-56.3±0.6 


14.1±0.9 


-49.3±0.5 






17.3±1.6 


34.9±0.3 


13±2 


37.6±0.6 


PSA 


5.54 


-65.54 


22.54 


-57.97 


-1.36 


-0.98 


22.30 


43.13 


10.08 


47.78 


AGS 




-65.97 


22.57 


-58.18 


-1.39 


-0.97 


22.49 


43.21 


10.13 


47.73 


HH 




-67.8±0.9 


21.3±0.7 


-58.6±0.3 






21.2±1.7 


45.2±0.5 


14±2 


51.5±0.5 


PSA 



TABLE VI: Same as Table IIV1 but for the low-fc potential. In this Table, only the AGS and HH results are reported. The 
phase-shifts and mixing parameters obtained by the recent PSA [35|] are also shown. 




FIG. 1: Differential cross section and neutron and triton analyzing powers A y o and Ao y for n- 3 H elastic scattering at E n = 1, 2, 
3.5, and 6 MeV neutron lab energies as functions of the cm. scattering angle. Results obtained using the AGS equation (solid 
lines), the HH expansion method (dashed lines), and the FY equations (dotted lines) using the I-N3LO potential are compared. 
For most of the cases the three curves coincide and cannot be distinguished. The experimental data are from Ref. [46( [. 
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FIG. 2: Same as Fig. [D but for p- 3 He elastic scattering at E p = 2.25, 4.05, and 5.54 MeV proton lab energies. The experimental 
data are from Refs. [29H33I . f35l] . 
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FIG. 3: Same as Fig. [Q but for the AV18 potential. 
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FIG. 4: Same as Fig. but for the AV18 potential. 



15 




FIG. 5: Same as Fig. [TJ but for the low-fc potential. Only the AGS and HH results are reported. 




FIG. 6: Same as Fig. [21 but for the low-fc potential. Only the AGS and HH results are reported. 
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FIG. 7: (Color online) Differential cross section, proton analyzing power, and 3 He analyzing power for p- 3 He elastic scattering 
at E p = 2.25, 4.05, and 5.54 MeV proton lab energies obtained using the I-N3LO (red bands), AV18 (blue bands), and the 
low-fc (cyan bands) potential models. The experimental data are from Refs. [29I4331 [35j | . 
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FIG. 8: (Color online) Same as Fig. [JJ but for the spin correlation A yy and A xx observables. The experimental data are from 
Refs. [US!. 




FIG. 9: (Color online) Same as Fig. [JJ but for the spin correlation A xz , A zx , and A zz observables (at E p — 5.54 MeV proton 
lab energy, only). The experimental data are from Ref. [3^ ]. 



